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TRANSVERSE  VIBRATIONS  OF  HOLLOW  THIN-WALLED 
CYLINDRICAL  BEAMS1 

By  Bernard  Bcdianskt  and  Edwin  T.  Kri  siewski 


SUMMARY 

The  variational  principle,  differential  equations,  and  bound¬ 
ary  conditions  considered  appropriate  to  the  analysis  of  trans¬ 
verse  vibrations  of  hollow  thin-walled  cylindrical  beams  are 
shown.  General  solutions  for  the  modes  and  frequencies  of 
cantilever  and  free-free  cylindrical  beams  of  arbitrary  cross  sec¬ 
tion  but  of  uniform  thickness  are  given.  The  combined  influence 
of  the  secondary  effects  of  transverse  shear  deformation,  shear  lag, 
and  longitudinal  inertia  is  shown  in  the  form  of  curves  for 
cylinders  of  rectangular  cross  section  and  uniform  thickness. 
The  contribution  of  each  of  the  secondary  effects  to  the  total 
reduction  in  the  actual  frequency  is  also  indicated. 

INTRODUCTION 

The  elementary  theory  of  bending  vibration  is  often  in¬ 
adequate  for  the  accurate  calculation  of  natural  inodes  and 
frequencies  of  hollow,  thin-walled  cylindrical  beams.  Such 
secondary  effects  as  transverse  shear  deformation,  shear  lag, 
and  longitudinal  inertia,  which  are  not  considered  in  the 
elementary  theory  of  lateral  oscillations,  can  have  appre¬ 
ciable  influence,  particularly  on  the  higher  modes  and 
frequencies  of  vibration.  The  effects  of  transverse  shear 
deformation  and  of  rotary  (rather  than  longitudinal)  inertia 
have  been  studied  by  many  on  the  basis  of  the  original  inves¬ 
tigations  of  Rayleigh  (ref.  1)  and  Timoshenko  (ref.  2). 
Anderson  and  Houbolt  (ref.  3)  have  presented  a  procedure 
for  including  the  effects  of  shear  lag  in  the  numerical  calcu¬ 
lation  of  modes  and  frequencies  of  box  beams  of  rectangular 
cross  section.  However,  there  does  not  appear  to  exist  a 
general  solution  for  the  vibration  of  hollow  beams  that  in¬ 
corporates  the  influence  of  all  the  secondary  effects  men¬ 
tioned. 

The  purpose  of  the  present  report  is  threefold:  First,  to 
exhibit  the  variational  principle,  differential  equations,  and 
boundary  conditions  appropriate  for  the  analysis  of  the 
uncoupled  bending  vibration  of  hollow  thin-walled  cylindrical 
beams;  second,  to  give  general  solutions  for  cantilever  and 
free-free  cylinders  of  arbitrary  cross  section  but  of  uniform 
thickness;  and  finally,  to  show  quantitatively  the  influence 


of  the  secondary  effects  by  means  of  numerical  results  for 
hollow  beams  of  rectangular  cross  section  of  various  lengths, 
widtlis,  and  depths. 
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SYMBOLS 

cross-sectional  area 
Fourier  coefficient 
effective  shear-carrying  area 
parameter  defined  in  equation  (30) 
constant 

modulus  of  elasticity 
shear  modulus  of  elasticity 
moment  of  inertia 

geometrical  parameter  defined  in  equation  (20) 
length  of  cantilever  beam,  half-length  of  free- 
free  beam 

parameter  defined  in  equation  (3K) 
maximum  kinetic  energy 
maximum  strain  energy 
half-depth  of  rectangular  l>eain 
half-width  of  rectangular  beam 
Fourier  series  coefficients 
integers 

frequency  coefficient, 

coefficient  of  shear  rigidity, 

coefficient  of  rotary  inertia,  J 

perimeter  of  cross  section 

distance  along  periphery  of  cross  section  (see 

fig-  1) 

wall  thickness 

longitudinal  displacement  in  z-direction 
vertical  displacement  in  g-diroction 
longitudinal  coordinate 
vertical  coordinate 

g -coordinate  of  center  of  gravity  of  cross  section 
shear  strain 
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«,  longitudinal  ■  train 

#  inclination  of  normal  with  vertical  (see  fig.  1) 

X  Lsgrangian  multiplier 

m  mam  of  beam  per  unit  length 

p  mam  density  of  beam 

•  longitudinal  direct  stress 

r  shear  stress 

u  natural  frequency  of  beam 

ut  natural  frequency  of  beam  calculated  from  ele¬ 

mentary  beam  theory 

tu  Kronecker  delta  (1  if  i=j;  0  if »« *j) 

9  constraining  relationship 

BASIC  EQUATIONS 

Assumptions.  The  problem  to  be  considered  is  that  of 
the  natural  bending  vibration  of  a  thin-walled  hollow 
cylindrical  beam  whose  cross  section  is  symmetrical  about  at 
least  one  axis  (see  fig.  1 ).  The  transverse  vibration  is  sup¬ 
posed  to  take  place  in  the  direction  of  this  axis  of  symmetry 
of  the  cross  section  so  that  no  torsional  oscillations  are 
induced. 

In  the  present  analysis,  the  following  simplifications  arc 
introduced: 

(a)  Changes  in  the  size  and  shape  of  the  cross  section  a  re 
neglected. 

(b)  Stress  and  strain  are  assumed  to  l>e  uniform  across  the 
wall  thickness. 

(c)  The  small  effect  of  circumferential  stress  upon  longi¬ 
tudinal  strain  is  neglected. 

In  accordance  with  statements  (a)  and  (b),  the  distortions 
of  the  vibrating  beam  are  completely  described  by  the 
vertical  displacement  w(x)  of  a  cross  section  and  the  longi- 


SymmetrcoJ 

obout 
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tudinal  displacement  u(x,»)  of  each  point  of  the  median  line 
of  the  beam  wall. 

The  longitudinal  and  shear  strains  are  given  in  terms  of 
u(r,«)  and  w(x)  as 


•‘“dr 


d«  .  dtr  .  „ 

y"  =  dH+d,  8,tt  9 

and  the  corresponding  stresses  Iterotnr 


where  8  is  the  inclination  of  the  normal  with  the  vertical  (see 
fig.  I  )• 

In  elementary  beam  theory,  when*  the  effects  of  all  shear 
distortion  are  neglected,  the  longitudinal  distortion  u(x,*)  is 
related  to  tile  vertical  displacement  trU)  In¬ 
i'  (s,*)^(y-y)  j" 

where  y  is  the  (/-coordinate  of  the  center  of  gravity  of  the 
cross  section.  In  the  present  report,  however.  uts.*)  is 
allowed  to  Ih*  perfectly  general,  so  that  shear  distortions  (and 
consequently  the  so-called  shear-lag  and  transverse-shear- 
deformation  effects)  are  fully  taken  into  account.  Further¬ 
more,  Ih-csusc  cross  sections  an-  not  constrained  to  remain 
plane,  the  inertia  effect  associate!  with  motion  in  the 
longitudinal  direction  is  mon>  properly  designated  as  the 
effect  of  longitudinal  inertia  than  the  effect  of  rotary  inertia. 

Variational  principle  and  geometrical  boundary  condi¬ 
tions. — The  variational  equation  to  be  written  is  appropriate 
to  beams  whose  ends  an*  either  fixed,  simply  supporter!,  or 
free.  For  some  such  beam  vibrating  in  a  natural  mode,  the 
maximum  strain  energy  is 

(5) 


»)  TypteM  eras  Met  ton.  (b)  Bi*n  conventions, 

U  QsaSgever  beam.  (d)  Free-free  beam 

VwillUB  ]. — Coordinate  systems  and  sign  conventions. 


where  u(x,»)  and  tr(/)  are  the  amplitudes  of  displaeement  for 
the  particular  mode  considered.  The  maximum  kinetic 
energy  is 

*+  j  f  pMv’d*  dx  (8) 

where  w  is  the  natural  frequency  of  the  mode  under  considera- 
tion  and  p  is  the  mass  density  of  the  beam.  The  second 
term  in  equation  (6)  constitutes  the  eontribution  of  longi¬ 
tudinal  inertia  to  the  kinetic  energy. 

A  natural  mode  of  vibration  must  satisfv  the  variational 
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equal  km 


t(U- n=o 


(7) 


Fixed  end: 


ir  0 


where  the  variation  is  taken  independently  with  respect  to 
*(/,*)  and  w(x)  and  witii  the  provision  that  both  «(/,«)  and 
ic(/)  must  satisfy  the  geometrical  boundary  conditions  of 
the  problem;  furthermore,  w(/,*)  must  be  periodic  in  the 
coordinate  *  with  a  period  equal  to  the  perimeter  p.  The 
geometrical  l«oundary  conditions  are  if =0  and  «=0  at  a 
fixed  end  and  only  tr=  0  at  a  simply  supported  end.  At  n 
free  end  no  geometrical  Imundarv  conditions  are  imposed. 

Differential  aquations  and  natural  boundary  conditions . 
Equations  (5).  («).  and  (7)  in  conjunction  with  the  usual 
procedure  of  the  calculus  of  variations  yield  the  following 
simultaneous  iutegrodifferential  equations  for  m  and  if: 

<>£[«<£+£■ »”•)>*<-'* -  ■» 


Simply  sup|M>rled  end: 


«  0 
if  » 


Free  end: 


E,  0 
0/ 


(|>  (it  (  sin  9^  sin  a  il*  II 


a:/*"  ■« 

d/ 


The  inlegriHlifrerenlial  equations  ih )  and  tit  >,  which  s|>ccify 
equilibrium  in  the  longitudinal  and  transverse  directions 
respectively,  can,  of  course,  be  written  directly  without 
recourse  to  the  variational  principle. 


where 


(J)  (it  (  ^  sill  'jfyj  s',r  *  )  d"  f  HU* tP  -  0 

p  -  pt  tlx 


(») 

(10) 


and  the  lannulary  equations  at  each  end  of  the  beam  are 

d>A.'f  tlx  -  0  (ID 

(it  (  ^  ^  sin  sin  6  tlx  it r=-  0  ( I  2) 


At  a  fixed  end.  Ixith  Imundary  equations  (II)  and  ( 12)  are 
satisfied  by  virtue  of  the  fact  that  the  geometrical  Ixiundary 
conditions  riapiire  that  Imth  in  and  ht  la*  zero.  At  a  simply 
supported  end  i>r—  I),  but,  since  Stiis.x)  is  perfectly  arbitrary, 
the  variational  process  forces  the  equality 

Et  ^"=0  (13) 

dr 

Finally,  at  a  free  end.  since  there  are  no  geometrical  con¬ 
straints.  both  in  and  iic  are  arbitrary  and  hence  the  varia¬ 
tional  proeess  forces,  in  addition  to  equation  (13),  the 
equality 

(jjf/f  (  ^  +  sin  9^  sin  9  tlx =  I)  (14) 

Equations  (13)  and  (14)  constitute  so-called  "natural 
lamiidary  conditions"  because  they  are  automatically  satis¬ 
fied  as  the  result  of  a  variational  process.  Equation  (13)  is 
recognized  as  the  eondition  of  zero  longitudinal  direct  stress 
while  equation  (14)  simply  stipulates  that  the  total  vertical 
shear  force  vanish. 

Thus  to  summarize,  the  appropriate  boundary  conditions 
required  for  the  solution  of  equations  (X)  and  (9)  are 


I  GENERAL  SOLUTIONS  FOR  CYLINDERS  OF  I’NIFORM  WALL 
THICKNESS 

The  following  exact  solutions  for  cylinders  of  uniform  wall 
thickness  arc  carried  out  by  means  of  Fourier  series  in  con- 
junction  with  the  application  of  the  variational  condition 
(eq.  (7)).  This  procedure,  which  does  not  require  explicit 
consideration  of  the  natural  lioundnry  conditions,  was 
believed  to  l»e  more  expedient  than  a  direct  attack  iq>on  the 
simultaneous  integrodilfcrential  equations  iX)  and  (9)  and 
;  all  their  associated  itoundury  conditions. 

Cantilever  beam.  The  geometrical  boundary  conditions, 
for  a  cantilever  beam,  as  previously  shown,  are 

if«))-  M  (!),*)  =  0 

j  (see  fig.  1).  Appropriate  assumptions  for  the  displacements 
j  ifU)  and  «(/.*)  are 

, ,  ,  «-■>  .  nr/ 

»'(/)  (  +  2L,  6,  COS  .  (lil) 

» . !  .:c:.  -  L 

j  and 

-L,  nir/  2  nr* 

«(/,*)  Zj  Zj  "■»>  s,n  .>i  ••ox  (Hi) 

m  rn  |„v.  H  p 

The  condition  «(()..<)  0  is  satisfied  by  each  term  of  equa¬ 

tion  ( Mi):  the  condition 

«•(»)  f+  v)  «  O’) 

I  »  .|_v. 

i 

is  introduced  into  the  variational  procedure  by  means  of  the 
Lagrangian  multiplier  method.  The  choice  of  the  imrticular 
trigonometric  functions  used  in  the  Fourier  series  (15)  and 
(It!)  was  guided  by  consideration  of  the  orthogonality 
j  required  for  the  simplification  of  expressions  in  the  strain 
j  energy.  The  constant  C  is  needed  in  the  expression  for 
V(/>  in  order  that  *r(/.)  be  unrestricted. 
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Using  equations  (15)  and  (16)  in  equations  (5)  and  (6) 
yields 


With  the  use  of  the  nondimensional  parameters 


r-r= 
1  /•£ 


1  fL(&  t?t(  m*  ">rx  2«r*Vj  .  , 

2  Jo  2ZCOS  2  L  C08  p  )d*d*  + 

1  ^  2»V  .  WITT  .  2»»s  . 

2  1  9  Gt  L£s  .S,'a-  T  *"*  2 L  “  >  + 

(.ii,6'co8  ^z+c)*,u  d'~ 

f  1,  t p/  cs„  •« sin  -*?  ^  2,rh dx  (,s) 

'l'o  make  equation  (18)  stationary  and  at  the  same  time 
satisfy  the  constrainin';  relationship 

*=^+^^.  =  0  (19) 

it  is  sufficient  to  sot 

«(r-7’-xv.)  =  o  (20) 

where  the  variation  is  with  respect  to  the  «’s,  ft's,  and  C 
considered  as  independent  variables;  here  X  is  a  Lagrungian 
multiplier.  This  variational  process  results  in  the  following 
equations; 

d(f  — 7)  i><p  ^  „  in*’  .1.  „(  iw\‘  A*L. 


k“  ~  ei  “ 

(26) 

t2-  E/ 

Ks  ~abgv 

(27) 

,  .  /  / 

-ptir'AV 

(28) 

JO-  16/ 

K  A*? 

(29) 

=«2-w(?y 

(30) 

equations  (21),  (22),  and  (28)  may  he  reduced  to 

inr*  Li  .  I  /iwV  I  .  I  ,  *, 

,5.,  4 ks*  As  A’  "  +  2  V27  it.,2  b,~‘2  b< 

i- i  2  L3\ 

((=1,3,5,  ..  .)  (31) 

(*s2fl,2+JC2/)(l+«0,)a<,+K*  P  A,ijb,-0 

(i  =  1.3, 5 _ ) 


(j=0, 1,2,  .  .  .)  (32) 

nt 

’+*■1-0  (33) 


~  n  dip  *  iiir’  .1.  „/ ixV  A,£ ,  For  j=0,  equation  (32)  becomes 

a*,  ~xaV  ,.4t.,6/  2  2-"‘’+n2Z)  2  6-~ 


=  0  (i=l,  3,  5,  .  .  .)  (21) 

d(Tj-7’)..  / (>  ,*£/>  r  /2jr\*Lp 

"  do,;  2L).  4  C+indo,,-*  6/^  J  )  4  «„  + 

'J'2 . 1  ;ft ,  -«V  1  +  i„()<7  „  =  0 

(i=l.  3.  5,  .  .  .) 

(j=0,  1.2,  .  .  .)  (22) 

W-T)  d*  -  2L.  trr,  % 

~a^ 


d,=-(^  sin  0 sin  ^ nwM  <i* 
pj  P 

As--  sin2  ft  <1* 


Equation  (34)  is  not  coupled  to  any  of  equations  (31)  to 
(33).  A  given  value  of  a,0  corresponds  to  the  amplitude  of 
the  ith  mode  of  longitudinal  oscillation,  and  if  this  value  of 
0(o  is  not  equal  to  0,  then  equation  (34)  simply  gives  the 
frequency  of  this  longitudinal  mode.  Consequently  those 
equations  in  equation  (32)  for  values  of  j=0  are  not  asso¬ 
ciated  with  transverse  bending  and  so  are  ignored  henceforth. 
For  the  remaining  values  of  j  (that  is,  jVO)  equation  (32) 
yields 

_~K'hA>li  .  (i=  1,3,5,...) 

”  kW+K'?  (j=  1,2,3,  .  .  .)  (35) 

Substituting  the  expression  for  a„  in  equation  (35)  into 
equation  (31)  and  solving  for  ft(  gives 


(-1  n  T3\ 

(-i  )-?k„’c+i* 


({=  1.3,5,  .  .  .  )  (36) 
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where 


„  tV  1  16  ka1  A,*4,1(tn)1  1 

8  lea*  7.  kJBS+IPn*  2 


In  the  appendix  this  expression  for  N,  is  shown  to  be  equivalent  to 


Nr- 


32 


B,*— 


tV  A 


16  -4S  2 


» V 


A.* 


(37) 


(38) 


Since  the  series  in  equation  (38)  is  considerably  more  quickly  convergent  than  that  in  equation  (37),  equation  (38)  should  be 
used  in  actual  numerical  calculations  of  Nt. 

Substitution  of  equation  (36)  into  equation  (33)  and  the  cons  training-relation  ship  equation  (19)  gives  the  following  two 
homogeneous  equations  in  C  and  X: 


A  M 

»r  JV.J 

o 

II 

►J  1*1 

(39a) 

© 

1! 

(39b) 

Finally  the  condition  for  a  nontrivial  solution  for  C  and  X  gives  the  frequency  equation 


iv  r i +*.*  ±  (~ y  4  ]  * +*v  ±  (-  n"*J  -  4 

L  * — Tj.s \nr/  .V, J  nr  N. 


2  1 


!+*.*  S  (-D  *  v- 

•  -1.3,5  Wf  A* 


.Sa  (V. 


=  0 


(40) 


which  the  frequency  parameter  ka  must  satisfy.  Since  the 
terms  of  the  infinite  series  which  appear  in  the  frequency 
equation  contain  kB  itself,  the  roots  of  equation  (40)  are  most 
conveniently  found  by  trial.  Fortunately  the  infinite  series 
in  equation  (40)  as  well  as  the  series  in  the  definition  of  Nt 
converge  rapidly  so  that  only  a  few  terms  are  needed  to 
evaluate  them  with  sufficient  accuracy. 

Once  kB  has  been  determined  for  a  particular  mode,  the 
corresponding  mode  shape  can  be  found  by  letting  C—  1  and 
solving  either  of  equations  (39)  for  X  and  then  finally  evaluat¬ 
ing  b,  and  a„  successively  from  equations  (36)  and  (35). 

Tree -free  beam — symmetrical  modes.— If  the  origin  of  a 
free-free  beam  of  length  2 L  is  taken  at  the  midspan  (see  fig. 
1 ),  the  form  of  the  Fourier  series  assumed  for  w(s)  and  u(i,») 
when  the  beam  is  undergoing  a  symmetrical  mode  of  vibra¬ 
tion  may  be  exactly  the  same  as  that  assumed  for  the 
cantilever  beam  of  length  L  (see  eqs.  (15)  and  (16)).  The 
only  difference  in  the  ensuing  calculations  is  that  the  con¬ 
straining  condition  (19)  is  not  introduced.  Consequently,  it 
can  be  readily  seen  that  the  frequency  equation  for  the 


symmetrically  vibrating  free-free  beam  is  obtained  from 
equation  (39a)  by  setting  X=0  and  is 


After  a  particular  root  kB  is  found  from  equation  (41),  the 
shape  of  the  corresponding  symmetrical  free-free  inode  may 
be  obtained  from  equations  (36)  (with  X  =  0)  ami  equations 

(35). 

Free-free  beam — antisymmetric*!  modea. — Consider  a 
free-free  beam  of  length  2 L  undergoing  antisymmctrical 
vibrations.  Explicit  consideration  need  be  given  only  to  the 
right  half  of  the  beam  (see  fig.  1 ).  and  for  this  half-beam  the 
only  geometrical  boundary  condition  that  must  be  imposed 
is  that  w(0)=0.  The  spanwise  displacement  u(0,»)  is  unre¬ 
strained  by  virtue  of  antisymmetry. 

Appropriate  assumptions  for  the  displacements  t»(x)  and 
u(r,»)  are  then 

ao  Ilf  r 

w(i)=^  ^  b.  sin  -  l  +Cx  (42) 
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and 
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_s,  mu  2  nr* 

U((,»)=  2^a,,.  cos  cos  (43) 

The  linear  portion  Cx  of  the  expression  forie(x)  is  needed  in  order  to  give  the  heam  sufficient  fmilom  at  the  lip  (x=  /„). 
The  choice  of  the  particular  trigonometric  function  in  the  series  expansion  for  «(x,*)  was,  as  in  the  east-  of  the  cantilever 
heam.  guided  by  consideration  of  the  orthogonality  required  for  the  simplification  of  the  expressions  in  the  strain  energy. 
'Hie  zeroeth  term  in  the  series  for  «(x,.«)  in  the  ((-direction  was  omitted  because  it  only  leads  to  the  frequency  equation  for 
longitudinal  oscillations. 

I'sing  equations  |42)  and  (43)  in  equations  (5)  and  (fit  yields 

i?  -  t  y  -  u.  &-■ r  -  rs  -  t+ 

sin  *(.£,  **  S  PWlT+r)J  h- 

p'“"Ls,,  c°s  "r  r°a  T)‘  ,44i 


The  variation  of  equation  (44)  with  respect  to  the  «'s. 
fc's,  and  (’  gives,  after  suitable  simplification. 

(Ki‘k*s+K*j*)a,l-K*fr  A,ijb,  ■-«  (»— 2.4.t>.  .  .  .) 

(j--=  1.2.3.  ..  .)  (45) 

£k-/*  -  kjkjk,,*  (  *  )'"]  an,—K‘  £  A,jf  (l 

(j=  1.2.3,  ...  I  (4ti) 

^  1  Lt  i H r‘  .  1  1  ,  1  .  ,i 

;i,  4  6-+2^A|- 


kH‘('L~ 0  (i  —  2,4,ti. 


. )  (47) 


a  l  Lt  .  <  L 

Ak*tA»"*AmUm~W 


kH‘  i:  -  (~l)""6„-f  IkSCL^tt  (48) 
.  "...  nr  i 

From  etpiation  (45) 

**  fL  A,ij 

«•>  ’  H  1  k  i  4-  fO  ?  b‘  •  -  •  >  (4!'> 

'  nj  ( 1 . 2. 3 _ I 

which,  except  for  sign,  is  the  same  expression  as  (hat  ob¬ 
tained  for  the  cantilever  and  symmetrically  vibrating  free- 
free  beams  (eq  (35)).  From  equation  (4ft) 

0n> “ HSka'+Ktj,V  ij~  l-2"’ - ]  <50) 

Substitution  of  equation  (40)  into  equation  (47)  gives 


*,=  -(-  I CL  (/=  2. 4.6....)  (51) 

i»  A i 


when-  A’i  is  defined  in  equation  (37). 

Substitution  of  equations  (.50)  and  (51 )  into  equation  l4S) 
and  simplification  gives  ms  the  frequency  equation  for  the 
antisymmetrically  vibrating  free-free  beam 

wDt O' iH *-*.4 **£*v+J]  • 


After  a  particular  value  of  kH  is  found  from  equation  (52 1, 
the  sha|M-  of  the  corresponding  ant isy linnet rical  free-free 
mode  may  la-  obtained  by  giving  C  the  arbitrary  value  of 
unity  and  calculating  the  b' s  and  rt’s  successively  from 
e<|iiations  (51),  (50),  and  (401. 

Discnssion  of  parameters.  The  parameters  entering  in 
the  frequency  equations  merit  discussion.  The  unknown 
natural  frequency  is  contained  only  in  the  frequency  co- 

/  fit 

efficient  l'n.  which  is  defined  by  the  formula  cc~  V  «> 

>  Mi¬ 
ami  is  in  eommoti  use  in  beam-vibration  analysis.  The 
parameters  ks  and  kK,  an-  identical  with  I  lie  shear  and 
inertiu  parameters  defined  in  reference  4.  which  considers 
the  effect  of  only  transverse  shear  and  rotary  inertia  on 
Ih-biii  vibrations.  The  quantity  .).«  which  np|M-ars  in  the 
present  definition  of  k„  is  actually  the  effective  shear-carrying 
area  when  plane  sections  are  constrained  to  remain  plane; 
that  is.  when  shear  lag  is  neglected.  The  remaining  param¬ 
eters  appearing  in  the  present  derivation,  namely.  .1  .ls, 
K,  and  .1|,  .12.  .  .  .  are  essentially  shape  parameters  which 
actually  depend  only  on  the  contour  of  the  cross  sect  ion ; 
as  shown  in  the  appendix. 

*V  =  .>  2  -V 

and 

•)  I  ^  I  > 

K‘  ,  V  52  ' 

»-l.*  .rni..i  i< 

and  the  .l.'s  are  simply  the  Fourier  coefficients  of  the 
function  sin  9.  which  is  dependent  only  on  the  shape  of  the 


I  • 

*■■■ 

i 
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cross  section.  These  shape  parameters  are  related  to  shear- 
lag  effects  and  their  interaction  with  transverse  shear 
and  longitudinal  inertia. 

The  effect  of  longitudinal  inertia  is  associated  » itli  the 
parameter  ku.  If  the  effect  of  longitudinal  inertia  is  to  l>e 
neglected,  it  is  sufficient  to  set  ku  equal  to  zero  in  the  final 
frequency  equation.  If  kKr  is  eq tial  to  zero,  H,  becomes  in¬ 
dependent  of  kt.  Appreciable  simplification  in  a  trial-and- 
error  solution  for  the  natural  frequency  then  results  since, 
with  /*,  independent  of  kB,  the  infinite  summation  contained 
in  .V,  is  also  independent  of  ka  and  need  he  calculated  only 
once  for  any  particular  beam.  As  is  shown  in  the  following 
section,  the  effect  of  disregarding  the  influence  of  longitudinal 
inertia  may  often  be  negligible. 

Without  presentation  of  details,  it  may  be  mentioned  that 
for  the  case  of  a  circular  cylinder,  which  has  no  shear  lag, 
all  the  .-l.'s  except  .h  vanish  ami  the  frequency  equations 
(40).  (41).  and  (.Vi)  may  be  put  into  closed  forms  identical 
to  those  given  in  reference  4.  Again,  if  in  the  general  fre¬ 
quency  equations  ks  is  set  equal  to  zero,  the  equations  may 
be  put  into  closed  forms  equivalent  to  those  of  reference  4 
where  only  rotary  inertia  is  considered. 

RESULTS  FOR  CYLINDRICAL  BEAMS  OF  RECTANGULAR 
CROSS  SECTION 

In  order  to  show  quantitatively  the  effects  of  shear  lag. 
transverse  shear  deformution,  and  longitudinal  inertia  on 
the  natural  frequencies  of  hollow  thin-walled  cylindrical  i 


beams,  numerical  calculations  have  been  performed  for 
cylinders  of  rectangular  cross  section  oscillating  as  free-free 
beams.  The  calculations  have  been  limited  to  symmetrical 
modes  of  vibration,  and  consequently  the  frequency  equation 
(41)  is  applicable.  For  ret  tangular  cross  sections  the  quan¬ 
tity  A',  may  be  put  into  closed  form  as  shown  in  the  appendix, 
ami  this  closed-form  version  of  A’,  was  used  in  the  calcula¬ 
tions.  A  value  of  K  (•'  equal  to  2.415  (appropriate  for 
aluminum  alloys)  was  assumed. 

The  results  of  these  calculations  art*  shown  in  figures  2,  .‘1, 
and  4.  In  figure  2.  the  ratio  of  the  nutural  frequency  u  to 
the  natural  frequency  obtained  from  elementary  beam 
theory  is  shown  as  a  function  of  the  plan-form  aspect  ratio 
Lb  for  cross-sectional  aspect  ratios  of  1.0.  il.ti,  ami  <*>. 
The  contribution  of  each  of  the  secondary  effects  to  the  total 
reduction  in  the  natural  frequency  for  the  cross-stud ional 

aspect  ratios  ^ --S.fi  ami  1.0  can  be  seen  in  figures  II  and  4, 

respectively  The  cross-sectional  aspect  ratio  of  ^  =  a* 

eorrosjtomls  to  the  limiting  case  of  a  beam  where  the  effects 
of  transverse  shear  deformation  and  longitudinal  inertia  are 
negligible  and  therefore  the  reduction  in  natural  frequency 

1  is  due  entirely  to  shear  lag. 

j  The  dashed  lines  in  figures  :i  and  4  show  the  reduction  in 
frequency  due  to  the  inclusion  of  the  effect  of  only  transverse 
shear  deformation  us  obtained  from  reference  4. 

The  long-  and  short-dash  lines  are  calculated  from  the 
frequency  equation  (41)  with  kHI  ()  ami  consequently 


(a)  First  symmetrical  mode.  (t>)  Second  symmetrical  mode.  (ci  Third  symmetrical  mode. 

Fiovrk  2. — Change  in  the  natural  frequency  of  a  symmetrically  vibrating  free-free  cylinder  due  to  the  inclusion  of  -ccoudary  effect*. 
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Mod*: 

Bt  symmetrical 


2d 


^  /('S'  ' M  symmetrical 


Effects  included: 
Transverse  shear  {ref.  4) 
Transverse  shear  and  sheor  lag 
Transverse  shear,  shear  log,  and 
longitudinal  inertia 


Bt  symmetrical 


2d  symmetricol  / 

y 


3d  symmetricol 


Effects  ^eluded 
Transverse  sheor  {ret  4) 
Tronsverse  shear  and  sheor  tog 
Transverse  sheor,  shear  log,  on<j 
longitudinal  mertio 


Kim'rf.  3  -  Contribution  of  transverse  shear  deformation,  shear  lag, 
and  longitudinal  inertia  to  the  reduction  in  natural  frequency  for 


Kiu  kk  4  Contribution  of  transverse  -diear  deformation,  -hear  lag. 
and  longitiniina)  inertia  to  tin*  reduction  in  natural  frequi-ncx  fur 

h  .  „ 


represent  (lie  reduction  in  natural  frequency  when  both. shear 
lag  and  transverse  shear  deformation  are  taken  into  account. 
Thus  the  hatched  area  between  the  dashed  and  flic  long-  and 
short-dash  lines  may  be  considered  as  showing;  the  additional 
reduction  in  natural  frequency  when  the  influence  of  shear 
lap  is  considered.  Finally,  the  solid  lines  are  calculated  with 
ku  taken  into  account,  and  consequently  the  shaded  area 
shows  the  additional  influence  of  longitudinal  inertia  in 
reducing  the  frequency. 

Examination  of  figures  3  and  4  and  the  etirves  for  ---  c°  in 
figure  2  shows  that  the  influence  of  shear  lag  increases  as  the 


cross-sect ionnl  aspect  ratio  increases:  whereas  the  influence 
of  transverse  shear  and  longitudinal  inertia  decreases  w  it li 
increasing  cross-sectional  aspect  ratio.  Indeed,  it  appears 

from  the  results  for  -  ~d.t>  llmt  for  this  aspect  ratio  the 
ii 

effects  or  longitudinal  inertia  may  already  be  considered 
practically  negligible. 

A  word  of  caution  concerning  the  interpretation  of  figures 
d  anti  4  may  be  in  order.  Since  in  some  eases  the  depth  of 
the  hatching  increases  w  it  It  increasing  /,  b.  it  might  appear,  at 
first  glance,  that  the  shear-lag  elfeet  increases  w  it  Ii  increasing 
plan-form  aspect  ratio.  However,  if  the  additional  effects 
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of  shear  lag  are  considered  on  a  percentage  basis  with  the 
dashed  line  as  a  base,  it  will  be  found  that  shear-lag  effects 
actually  reduce  in  percentage  with  increasing  Ljb.  A  similar 
criterion  should  be  used  in  judging  the  influence  of  longi¬ 
tudinal  inertia. 

CONCLUDING  REMARKS 

The  numerical  calculations  show  that  secondary  effects 
have  appreciable  influence  on  the  natural  frequencies  of 
rectangular  box  beams  of  uniform  wall  thickness.  These 
results  constitute  an  indication  of  the  probable  inadequacy  of 
elementary  beam  theory  for  the  vibration  analysis  of  actual 
aircraft  structures  of  the  monocoquc  and  semimonoeoqile 


type  and  emphasise  the  need  for  practical  calculation  pro¬ 
cedures  for  such  structure*  that  would  take  into  account 
transverse  shear  deformation,  shear  lag,  an«l,  w hen  necessary, 
longitudinal  inertia.  The  general  solutions  presented  for 
cylinders  of  uniform  thickness,  as  w  ell  as  the  numerical  results 
for  rectangular  liox  I  teams,  should  Ite  useful  in  the  assessment 
of  the  accuracy  of  any  procedure  of  this  kind  that  may  Ite 
developed. 

laxoLEv  Akko.yai'tk'ai.  Lai«>kat<>k>, 

National  Advisory  ('omuittll  for  Aliionai  tics. 

Lanolkv  Field,  Va.,  January  Jl,  19M. 


APPENDIX 


TRANSFORMATION 

Expressions  for  As/A,  I,  and  K2. — If  sin  0  is  expanded  into 
a  Fourier  series 


.  ■  2/its 

sm  0=  Tj  An  sin 

n*l,5,3  p 


(Al) 


the  Fourier  coefficients  /l,  are  the  same  os  those  defined  in 
equation  (24) ;  that  is, 


2  /  ■  .  .  2/tas  . 

,=  -(17  sm  0  sin  (Is 
PJ  P 


(A  2) 


The  effective  shear  area  A3  (eq.  (25))  can  now  be  written 
as  a  function  of  the  Fourier  series  expansion  for  sin  0  as 

t  (  £  ^  A.  sin  (A3) 

With  the  use  of  the  appropriate  orthogonality  conditions, 
equation  (A3)  becomes,  after  the  integration  is  performed, 

Aa—~  2  -V=4  £  .4.* 

2  Rrnb  2  n,  1,2.3 


or 


_ 1  v  A  2 

A  t)  ZmJ  ^ " 


(A4) 


The  moment  of  inertia  /  of  a  cylinder  is  defined  as  (see 

fig-  1) 

7=J*  yHds—Ay1  (A5) 

where  y  is  the  j/-distance  to  the  center  of  gravity  of  the  cross 
section  and  is  given  by 

r  yt  <h 
-  .ft) 


v- 


But 


pt 

y—J  sin  6ds 


(A6) 

(A7) 


of  parameters 

or 


£  -«./  (i-co*2""*)  (as, 

2m  x  \  p  / 

and,  consequently, 

2ntr 


(A9) 


With  the  use  of  equations  (AS)  and  (All),  tile  expression 
for  /  in  equation  (A5)  becomes 


r  -I  P  *t~\  d, 
id 


(A  10) 


With  the  series  expansion  for  /  in  equation  (Alt)),  the  param¬ 
eter  K1,  as  defined  in  equation  (29),  becomes 


IT  a  I  n  n*  1,2,3  W 


(All) 


Transformation  of  expression  for  N,. — In  equation  (37) 
N,  was  defined  as 

r_iV  iV  ™  A  ^  id  A,*  1  ,  ,  . 

4'  ‘~~8ks*  IMs1  Asn&,tka1Bl'+Kini  i  "  ^  ' 


The  infinite  series  that  appears  in  this  expression  converges 
as  AJ  and  therefore  is  a  relatively  slowly  converging  series. 
In  order  to  increase  its  rate  of  convergence,  the  following 
transformations  arc  made. 

By  adding  and  subtracting  A J/K2  inside  tire  infinite 
summation  in  equation  (Al2)  and  using  equation  (A4),  the 
equation  simplifies  to 

As  WftW-i  »**  (A,3) 

By  adding  and  subtracting  AJIK‘n‘  inside  the  infinite 
summation  in  equation  (A  13)  and  using  equation  (All),  the 
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expression  for  -V,  can  be  transformed  to 


-V, 


32 


Bf- 


16  A,  ‘ ts  .h.tK'n>(ks’Hl'Ak'n*)- 


(A14) 


The  infinite  series  in  equation  (A  14)  converges  as  A.’/n4  and  therefore  is  considerably  more  quickly  convergent  than  the 
aeries  in  equations  (Al2)  and  (Al3),  which  converge  as  A.1  and  A, 'In’,  respectively. 

Closed  fora  of  N,  for  cylindrical  beams  of  rectangular  cross  section.  -For  a  cylindrical  beam  of  rectangular  cross 
section,  with  dimensions  as  shown  in  figure  2,  it  is  possible  to  write  the  expression  for  X,  in  a  closed  form.  The  param¬ 
eters  for  such  a  cross  section  become 

A,=4«f  A 

A  =  4  (a  +  6)1=  pi  I 


A.  =  0 


4  2  nwb 

- -  C08 

nr  p 


(n  even)  " 
(n  odd) 


(A  15) 


With  equations  (A15)  the  parameter  A',  shown  in  equation  (A  12)  becomes 


or 


,  2nwb 

v  _  lV  _  *  P  P 

*'  s ks’  4 kfamfrUk.\ut  , 


•>  r» 


(4nr4  \ 

V'  *  ,  .  C°*  P  I  I  ,  , 


(Aid) 


(A  17) 


Each  of  the  infinite  summations  in  equation  (Al7)  ran  now  be  written  in  closed  form  as  shown  in  reference  5,  and  the 
closed  expression  for  A’<  then  becomes 

-i*.1  (A18) 
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